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ABSTRACT 

Young massive stars in the central parsec of our Galaxy are best explained by star 
formation within at least one, and possibly two, massive self-gravitating gaseous discs. 
With help of numerical simulations, we here consider whether the observed popula- 
tion of young stars could have originated from a large angle collision of two massive 
gaseous clouds at i? ~ 1 pc from Sgr A*. In all the simulations performed, the post- 
collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric 
outer filaments, consistent with the observations. Furthermore, the radial stellar mass 
distribution is always very steep, oc again consistent with the observations. 

All of our simulations produce discs that are warped by between 30° to 60°, in ac- 
cordance with the most recent observations. The 3D velocity structure of the stellar 
distribution is sensitive to initial conditions (e.g., the impact parameter of the clouds) 
and gas cooling details. For example, the runs in which the inner disc is fed inter- 
mittently with material possessing fluctuating angular momentum result in multiple 
stellar discs with different orbital orientations, contradicting the observed data. In all 
the cases the amount of gas accreted by our inner boundary condition is large, enough 
to allow Sgr A* to radiate near its Eddington limit over ^ 10^ years. This suggests 
that a refined model would have physically larger clouds (or a cloud and a disc such 
as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 parsec 
as in our simulations. 

Key words: accretion, accretion discs - Galaxy: centre - stars: formation - galaxies: 
active 



1 INTRODUCTION B-type stars l|Ghez et al.. 200dGhez et al.. 20051 ). These 

stars (commonly referred to as S-stars in the Galactic 
The power output of the central parsec of our Galaxy r-. ^ n j_ i j_i .l • i 

f ' Centre community) appear to be on rather eccentric orbits 

is dominated by young (aged 6 ± 1 million years i- ^ m ^ i -^i ^ ^ ^ i i ^ i 

~^ ' " — ; 1 , , p distributed without any notable preference to a plana r 

massive He-r stars flKrabbe et al.. 199a ). Many ol c ^- / \t^- i r~n — onndr 7 — \ — nnnA 

,^ \ ^ \ , , configuration (e.g., lEisenhauer et al.. 200aLu et al.. 20061 ). 

these stars are located m a well-denned thm disc . , ^ ■ ^ i -r o ^ i • n 

At present it is not clear it the b-stars are physically 

ated to the more massive "disc" -stars, as the initial 

3its of the S-stars could have been very different from 

their present day ones due to resonant relaxation processes 

ond more diffuse disc rotating counter clockwise ttz j — rr^^ ; r,r,r,A rm • ii • in 

r- ; ; — „ |, ; ; — " i, , , , flfiopman fc Alexander. 200q . thus, in this paper we shall 

jGenzel et al.. 2003Paumard et al.. 20061 ) although ^ ' , , — ^"j — f{, ' + j 



structure that rotates clockwise as seen on the sky ^ 

l|Levin fc Beloborodov. 200jGenzel et al.. 200dPaumard et al.. 2ob'^^ .^^ 
The remaining stars can be classified as a sec- 



concentrate on the puzzle of the young massive stars found 
the statistical significance of this feature is disputed ^ • i • i j_i ^ -^.i 

r; ; -i ° outside the inner arcsecond as these present us with a 

IILu et al.. 20061 ) . Stars appear to be on more eccentric orbits , i , , 



m 



— ; '"^ ~'' T ^ cleaner laboratory environment. 

the second disc (Paumard ct al., 200c Lu et a l.. 20061 ). 



The bright He-I stars seem to be excluded from the central iGerhard (200ll ) suggested that the stars may 

arcsecond (1 ^ 0.04 pc when viewed from the « 8 kpc ^^^^ ^^^^ ^^^^^^ ^ ^j^^^^^^ ^^^^ p^^^^^^ 
distance to the Galactic Centre (GC)), which instead ^ ^^^.^^ ^^^^ ^1^^^^^ ^j^^ ^1^^^^^,^ ^^^.^ ^^^1^ 

contains at least a dozen less massive but still quite young ^^^^ ^^^^^ ^^^^^^-^ dynamical friction with the back- 

ground stars, and eventually the stars would form 

a disc populated by stars on eccentric orbits (e.g., 

* E-mail: alexander.hobbsaastro.le.ac.uk [Hansen fc Milosavlievic. 200dKim et al.. 2004Gurkan fc Rasio. 200dLevin et 
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However, infrared observation s find no young massiv e stars 
outside the central ~ 0.5 pc (jPaumard et al.. 20061 ) which 
should be present in this model. In a ddition, Chandra X-ray 
limits ()Navakshin fc Sunvaev. 2005^ constrain the number 
of X-ray emitting Young Stellar Objects to well below what 
is expected in the cluster disruption scenario. 

An alternative model is the in-situ model in which 
stars form inside a s elf-gravitating gaseou s disc (e.g.. 



happened, e.g. a Giant Molecular Cloud (CMC) with a sub- 
parsec impact parameter (in relation to Sgr A*) could have 
self-colli ded and become partially bo und to the central par- 
sec (e.g.. lNavakshin fc Cuadra. 20051 '). Alternatively, a GMC 
could have struck the Circumnuclear Disc (CND) located a 
few parsecs away from Sgr A* , and then created gas streams 
that settled into the central parsec. 

It is very hard to estimate the probability of a 



iLevin fc Beloborodov. 2003JSIavakshin fc Cuadra. 200EtPaumard etlaild-3fl0fllfl collision in the central part of the Galaxy. 
In fact, theorists had expected gaseous massive discs around iHasegawa et al. (igQ^D estimated the rate of GMC collisions 

supermassive black holes to form stars or plane ts (e.g., as one per 20 Myrs in the central 100 pc, but the esti- 

iPaczvriski. 197dKolvkhalov fc Sunvaev. 198 C Shlos man fc Begelma na-ttSSS CDlltat&riZhbre,tlfiSS GammidcfiOffl fflafiditiidbq M)0^ 



long before the properties o f the young stars in t he GC 
became known. Recently, iNavakshin et al. (2007^ have 
numerically simulated the fragmentation process of a 
geometrically thin gaseous disc of mass > 10" Mq for 
conditions appropriate for our GC (albeit with a rather 
simple cooling prescription) and f ound a top-heavy mass 
function for the stars formed there. [Alexander et al. (20081 ) 
extended numerical studies to the fragmentation of eccentric 
accretion disks. 

However, the in-situ model has not so far addressed in 
detail the origin of the gaseous discs themselves, although 
similar discs are believed to exist in AGN and quasars. In- 
deed, sub-parsec massive gaseous discs are invoked as means 
of feeding superma ssive black ho les' immense appetite for 
gaseous fuel (e.g. iFrank et al.. 2002 ). Such discs are actu- 
ally observed as the site s of maser emission i n some of the 
nearby galactic centres l|Mivoshi et al.. 1995h . But Sgr A* 
is currently not a member of the AGN club, let alone 
the even more powerful and elitist quasar community. In- 
deed, Sgr A* bolometr ic luminosity is around ^ 10~^ of it s 
Eddington limit (e.g.. iNaravan. 200jBaganoff et al.. 2003t ). 
The amount of ionised gas currently present in the inner 
arsec is estimated at perhaps a mere hundred solar masses 
Paumard et al.. 20041). There are also tight observat ional 
constraints (|Falcke fc Melia. 199'iCuadra et al.. 2003h on 
the presence of an optically thin disc on smaller scales (say 
0.01 pc or less). 

We believe that the requisite gaseous discs could not 
have been assembled by viscous transport of angular mo- 
mentum as is expected to be the case in other systems 
l|Shakura fc Sunvaev. 197jFrank et al.. 20021 ). The viscous 
time scale of a thin, marginally self-gravitating disc can be 
estimated as tvisc = ( Mbh/Mdisc)^(af^)~^ , where Mbh ~ 
3.5 X lO'' M0 and Afdisc are masses of the blackhole and 
the disc, respectively, and Q is the Keplerian angular fre- 
quency. For a = 0.1 and Mdiac ~ lO" Mp (reasons f or 
this choice of mass are discussed in lNavakshin et al. (20061 ) ). 
tvisc ~ 3 X 10^ years at r = 0.0 3 parsec, the inner edge of 
the discs l|Paumard et al.. 20061 ). and 10^ years at 0.3 par- 
sec. These times are very long compared with the age of 
the stellar systems. In fact, the gaseous discs would have 
to evolve even faster as t he two stellar system s are co-eval 
within one million years (jPaumard et al.. 2006h . 

It is also not clear that one should expect any sort of a 
viscous quasi-steady state for the gas in our Galactic Centre. 
There is no strong evidence for other similar star formation 
events in the central parsec within the last ~ 10* years. 
Therefore, the one-off star formation event appears to be 
best explained by a one-off deposition of gas within the cen- 
tral parsec. There are several ways in which this could have 



ther direction due to the lack of knowledge about the size 
distribution of GMCs. However, there is observational evi- 
dence that GMCs can be put on orbits significantly differ- 
ent from the sim ple circular orbits seen in the inner Galaxy. 
IStolte et al. (2008) recently measured the proper motions of 
the Arches star cluster. This star cluster is about 2.5 mil- 
lion years o ld and has a mass in excess of lO" M© (e.g., 
iFiger. 20041 ). It is presently about 30 pc away from Sgr A* 
in projection. Its orbit must be strongly non-circular as its 
3D velocity is over 200 km/sec in the region where the cir- 
cular velocity is only ~ 110 km/sec. Therefore, this suggests 
that massive GMCs can be on highly non-circular orbits. It 
does not seem implausible that one of these clouds would 
pass within a few parsec of Sgr A* . 

In this paper we explore such a one-off collision event 
in a very simple setup O. We allow two massive, uniform 
and spherical clouds on significantly different orbits to col- 
lide with each other one parsec away from Sgr A*. The re- 
sulting gas dynamics, in particular the way in which gas 
settles into the inner parsec, is the focus of our effort here. 
We find that the collision forms streams of gas with vary- 
ing angular momentum, both in magnitude and direction. 
Parts of these streams collide and coalesce to form a disc, 
and remaining parts form one or more orbiting filaments. 
As the gas cools, it becomes self-gravitating and stars are 
born, usually in both the disc and the filaments. This overall 
picture is discussed in ij4.1l A reader mainly interested in a 
comparison between the simulations and observations may 
find sect ions [5] [6] and [7] most relevanlQ. In ^7.2\ we argue that 
our simulation results suggest that the location of the orig- 
inal impact was somewhat further away from Sgr A* than 
we assumed here, e.g., perhaps a few parsecs. 



2 NUMERICAL METHODS 

The numerical approach an d the code u sed in this paper 
is the same as in .Nayakshi n et al. (20071 ) with only slight 
modifications. We employ GADGET-2, a smoothed par - 
ticle hydrodynamic (SPH)/N-body code (|Springel. 20051 ). 
The Newtonian N-body gravitational interactions of parti- 
cles in the code are calculated via a tree algorithm. Artificial 
viscosity is used to resolve shocks in the gas. Gas cools ac- 
cording to du/dt = —u/tcooi{R), where cooling time depends 
on radius as 



(1) 



^ For movies of the simulations the reader is directed to 
|http : //www, astro . le ■ ac .uk/~aphll/movles .htmll 
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where tdyn = and = {GM\,\^/ R^Y^'^ is the Keplerian 
angular velocity, and /3 is a dimensionless parameter. This 
approach is motivated by simulations of marginally stable 
self-gravitating gas eous discs, where /3 is expected to be of 
the order of a few l|Gammie. 200]|Rice et al.. 2005h . In the 
problem at hand, 13 would most likely be much smaller than 
this during the collision phase of the clouds as gas can heat 
up to high temperatures (see Sj7]2}, but it could be much 
larger than unity when gas cools into a geometrically thin 
disc or a thin filament. In the framework of our simplified 
approach, and given numerical limitations, we only consider 
two values for 13. The (3 = 1 tests approximate the situa- 
tion when cooling is inefficient, whereas [3 = 0.3 cases corre- 
spond to rapid cooling. Each of our runs used between 3 to 
6 months of integration time on 12 processors in parallel. 

Gas moves in the gravitational potential of Sgr A*, mod- 
elled as a motionl ess point mass with Mbh = 3.5 X 10^ Mq 
l|Reid et al.. 19991 ) at the origin of the coordinate system, 
and a much older relaxed isotropic stellar cusp. For the lat- 
ter component we use the stellar de nsity profile derived f rom 
near-IR adaptive optics imaging bv lGenzel et al. (20031) : 




R< R cusp 

R > i?cusp 



(2) 



where 7?cusp = 10 arcseconds 



0.4 pc), and p^o = 1.2 x 

10** Mg/pc^ 

Note that the enclosed stellar cusp mass is ap- 



proximately 6 X 10 M0 at R 



Rc 



and equals 



Sgr A* mass at i? ~ 1.6 pc. We do not include 
the gravitational potential of the putative stel- 
lar mass black hole cluster around Sgr A* (e.g., 
iMorris. 199 jMiralda-Escude fc Gould. 2000l ). as at present 
its existence and exact properties are highly uncertain 
l|Freitag et al.. 200dSchodel et al.. 200iDeegan fc Navakshin. 

Star formation and accretion in SPH is tradi- 
tionally modelled via the "sink particle" formalism 
l|Bate et al.. 19951 ). where gravitationally contracted gas 
haloes are replaced by the coUisionless sink particles of same 
mass. We follow the same recipe, allowing Sgr A* to accrete 
all gas particles that are within -Race ~ 0.1 arcsecond (ex- 
act values vary between different cooling parameters, cf. f|3]). 
This is significantly smaller than the inner radius of the ob- 
served "disc" populati on of young massive stars, i.e., ~ 1" 
jPaumard et al.. 20061 ). 

To model star formation, new sink particles are intro- 
duced when gas density exceeds 



Pcrit = PO + j4co1PBH 



(3) 



where po = 1-67 x 10 g cm ^, pbh is the tidal den- 
sity, an d Acoi = 50 (this method is based on that of 
[Nava ks hin et al. (200i ) , which found that the results are 
not sensitive to the exact values of po and A^d, provided 
they are sufficiently la rge). For reasons fully explained in 
iNavakshin et al. ^2007^ , we distinguish between "first cores" 
- star particles with M < Mcoic = 0.1 Mq, and "stars" with 
M > Mcoro. Physically, first cores represent massive gaseous 
clumps that have not yet collapsed to stellax densities, and 
thus have finite sizes that are set to 7?core = lO^** cm. First 
cores are allowed to merge with each other if they pass within 
a distance d < 2i?corc. 



In contrast, stars are not allowed to merge, but can ac- 
crete first cores if they pass within a distance of Rcora- Ac- 
cretion of gas onto stars and proto-stars is calculated using 
the Bondi-Hoyle formalism, with the accretion rate given by 

(GMf 



M = 47rp 



(Av2 + c2)3/2 



(4) 



where M is the mass of the sink particle, p is the ambient gas 
density, Cs is the sound speed, and Av is the relative velocity 
between the gas and the accreting particle. The accretion 
rate is capped at the Eddington accretion rate: 

inrUpR^c 



MEdd = 



(5) 



where nip is the proton mass, c the speed of light and ctt 
the Thomson scattering cross-section. Note that this ex- 
pression depends only on the size of the object, i?, . For 
first cores, _R, — -Rcoro, which yields a very high accre- 
tion rate limit of almost a solar mass per year for 7?core = 
10^** cm. For stars, we use the observational results of 
iDemircan fc Kahraman" fl99]|Gorda' fc Svechnikov (19981 ): 



1.09 



1.29 



M 
M 



for M < 1.52 M, 



©1 



for M > 1.52 M, 



(6) 



(7) 



R 

Rq 

R 

Rq 

For R — Rq , this yields an Eddington accretion rate limit of 
~ 5 X 10~* Mq year"'^ . The actual gas particles that are then 
accreted by the sink particle are chosen from its neighbours 
via the stochastic SPH method o f lSpringel (20051 . The mass 
and momentum of the SPH particles accreted are added to 
that of the respective sink particle. 

Finally, we shall acknowledge that this plausible star 
formation prescription cannot possibly do justice to the full 
complexity of the physical star formation process. There- 



2007D, re, the stellar mass function obtained in our simulations 
should not be trusted in detail. However, our focus is on 
the dynamics of gas as it settles into a disc or a filamentary 
structure. When star formation does occur, the gas tem- 
perature is very much smaller than the virial temperature 
in the simulated region ( ^ 10^ A"). Gas orbits thus deter- 
mine the resulting stellar orbits which we compare to the 
observed ones. We believe that this important aspect of our 
simulations is reliably modelled. 

The set of units used in th e simulations is Mu = 3.5 x 
lO'^ Mq, the mass of Sgr A*(e.g.. [Sch"odel et al.. 2003 ). i?„ = 
1.2 X 10^*^ cm ~ 0.04 pc, equal to 1" when viewed from the ~ 
8 kpc distance to the GC, and tu = l/fl{Ru), the dynamical 
time evaluated at Ru, approximately 60 years. We use R 
to signify distance in physical units and the dimensionless 
r = R/ Ru throughout the paper interexchangeably. 



3 INITIAL CONDITIONS 

Simulations were performed for several sets of initial con- 
ditions, all of which comprise a collision between two gas 
clouds at the edge of the inner parsec of the GC. 2,625,000 
SPH particles were used in each simulation. The specific pa- 
rameters for each run, labeled SI to S6, can be found in 
Table 1. Each cloud is spherically symmetric and of uniform 
density, containing less than 1% of the mass of Sgr A* . The 
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clouds are composed of molecular hydrogen, with the mean 
molecular weight set to jj. — 2.46. The initial temperature of 
the clouds is set to 20K. A more complex model could have 
included a turbulent velocity and density field, but for prac- 
tical purposes we limit ourselves to a smaller set of input 
parameters for this first study. 

We define the primary to be the larger cloud, with a 
radius of Ri = 0.2 pc and Mi = 3.4 x IO^Mq, whilst 



the secondary has a radius of R2 — 0.172 pc with A/2 = 
2.6 X 10* M0. The initial positions of the cloud centres are 
the same for all the simulations, and are given in Table 1 
in dimensionless units. The initial density of the clouds is 
slightly above the tidal shear density at their initial loca- 
tions, so the clouds are marginally stable to tidal shear by 
Sgr A*. 

In all the tests the initial velocity of the primary cloud is 
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kept fixed, and is of tfie order of tfie Keplerian circular veloc- 
ity for that radius. Orbits in a cusped potential are however 
not Keplerian. The initial specific energy and angular mo- 
mentum I of a particle can be used to find the pericenter, rpc, 
and the apocenter, rap, of the orbit. The orbit's eccentricity 
e is then defined via 



In our simulations particles occasionally obtain hyperbolic 
orbits. For analysis (plotting) purposes only, the eccentrici- 
ties are capped at 1. In the stellar potential used here (equa- 
tion [2|, the orbit of the primary cloud is slightly eccentric, 
with pericentre and apocentre of 25 and 31.5 respectively, 
and with eccentricity e = 0.12. 

The initial trajectory of the secondary cloud is varied 
between the tests to cover a small range of possibilities. The 
parameters for this trajectory in terms of pericentre, apoc- 
entre and eccentricity are given in Table [T] The collision 
itself is highly supersonic, as is expected for molecular gas 
falling into the potential well of a supermassive blackhole. 

We ran tests with cooling parameter l3 — 1 and /3 = 0.3 
(see equation [TJ. These values are low enough so that frag- 
mentation would occur if and when regions of the ga s 
became self-gravitating l|Gammie. 200]|Rice et al.. 20051 ). 
Since the faster cooling runs were expected to require on 
average shorter timesteps the accretion radius for (3 — 0.3 
was set to race = 0.33 whilst for /3 = 1 a smaller value of 
race = 0.06 was used. 

We have performed convergence tests (at the referee's 
request) on the numerical viscosity parameter a and on 
Naph- We find that there is good convergence on the dis- 
tribution of gas when varying either of these parameters, 
although there are minor differences within the innermost 
part of the computational domain, namely the inner 0.5"-l" 
around the central black hole. These differences are small 
and can be explained by a different amount of gas being ac- 
creted onto the black hole when either viscosity or resolution 
are changed. 



4 SIMULATIONS 
4.1 Overall picture 

In all of our simulations, the clouds undergo an off-centre 
collision at time f ~ 10 in code units (~ 600 years). As 
the cooling time is longer than the collision time, fcoii ~ 
{Ri + iZ2)/(|vi — V2I), the clouds heat up significantly and 
hence initially expand considerably. This thermal expansion 
also modifies the velocity of the different parts of the clouds 
by giving the gas velocity "kicks" . The net result is a distri- 
bution of gas velocities that is much broader than what one 
would get if the two clouds simply stuck together, i.e. in an 
inelastic collision. 

The collision and resultant mixing of the clouds leads 
to angular momentum cancellation in some parts of the gas. 
Regions of gas that acquire smaller angular momenta infall 
to the respective circularisation radius on the local dynam- 
ical time. A small scale disc around the black hole is thus 
formed on this timescale. Regions of the clouds that did not 
directly participate in the collision are affected less and re- 
tain more of their initial angular momentum. These regions 



are initially not self-gravitating, and so are sheared by the 
tidal field of Sgr A* . They later cool and result in filaments 
of length comparable to the initial sizes of the clouds' orbits. 

As time progresses, parts of the filaments collide with 
each other or with the inner disc if the pericentres of their 
orbits are small enough. The inner disc therefore gains mass 
in an asymmetric non-planar manner, leading to the disc 
becoming warped and changing its orientation with time. 

Eventually, due to gravitational collapse, high density 
gas clumps are formed. Sink particles are introduced inside 
these clumps and are allowed to grow in mass via gas ac- 
cretion and mergers with first cores as explained in !j2] Star 
formation takes place in all of our simulations. Although the 
precise locations and the stellar mass distributions differ we 
always end up with stellar populations forming in two dis- 
tinct locations - in the disc around the black hole and in the 
outer filament (s). Disc populations generally have top-heavy 
mass functions whilst filament populations always consist 
entirely of low-mass stars; however, as noted in 321 the mass 
spectrum of st ars formed depends on a very simple prescrip- 
tion (e.g., see iNavakshin et al.. 2007h . We expect the mass 
spectrum of our models to change when more realistic cool- 
ing and feedback physics is included. 

4.2 SIMULATION SI 

4-. 2.1 Gas dynamics 

Figure [1] presents four snapshots of simulation SI showing 
gas column density and positions of the stars viewed along 
the z-axis. 

The initial mass deposition into the central few arcsec- 
onds forms a small-scale disc that changes its orientation 
by about 30° with time, stabilising by t ~ 300. As a result 
of the small impact parameter of the cloud collision, the 
secondary cloud is largely destroyed and does not collapse 
into a coherent filament. Instead, the diffuse gas left over 
from the impact gradually accretes onto the disc that forms 
around the black hole, growing it in mass and radius. The 
inner disc therefore ends up extending out to r ~ 10 in code 
units. The gas left over from the primary cloud collapses into 
a filament which orbits at a larger radius, separate from the 
disc. 

It is useful to analyse the re-distribution of the gas as a 
result of the collision in terms of orbital parameters, namely, 
angular momentum and eccentricity. Gas with specific angu- 
lar momentum Z = |r x v| circularises at the circularisation 
radius, rc, given by the relation 

P = GM(rc)rc , (9) 

where M(r) is the total enclosed mass within radius r. Fig- 
ure [2] shows the distribution of gas particles as a function of 
the circularisation radius of their orbits as a histogram for 
several snapshots. At time t — 0, the distribution is highly 
peaked around the clouds' initial circularisation radii. Due 
to the collision, however, a tail to small rc appears in the 
distribution very quickly. This is not due to viscous angular 
momentum transport, as in thin accretion discs, but is rather 
due to cancellation of oppositely directed angular momen- 
tum components in shocks. Some of the gas particles aquire 
particularly small angular momenta and hence fall to small 
radial separations from Sgr A* on nearly radial trajectories. 
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Table 1. Initial conditions of the simulations presented in the paper. The meaning of the symbols in the Table are: /3 is the cooling 
parameter, ri, r2, vi, V2 are the initial positions and velocity vectors of the two clouds, respectively; rpe and rap are the pericentres 
and the apocentres of the two clouds; e is their orbit's eccentricity; 6 is the angle between the orbital planes of the clouds, and b is the 
impact parameter. 
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/3 


ri 


r2 


Vl 


V2 


|vi - V2 
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29.1 


0.12 


0.38 


151 
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(0,-0.11,-0.21) 
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(22,6,7) 


(0,0.2,0) 


(0.16,-0.11,-0.21) 


0.41 


25 


31.5 


21.8 


43.1 


0.12 


0.33 
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2.4 


S6 


1 


(25,0,0) 


(22,6,7) 


(0,0.2,0) 


(0.16,-0.21,-0.11) 


0.45 


25 


31.5 


21.5 


43.4 


0.12 


0.34 


147 


5.3 




[0.04pc] r[arcsec] 



Figure 2. The distribution of SPH particles over the circularisa- 
tion radius of their orbits. The black, cyan, blue and red curves 
correspond to times t = 0, 25, 50, 500, respectively. As the result of 
angular momentum cancellation due to shocks, particles "leak" to 
small radii, establishing there a disc that eventually forms stars. 



Figure 3. Gas eccentricity defined on radial shells as a func- 
tion of the shell's radius at different times. Black, blue and red 
curves correspond to time t = 50, 100 and 250. Note that initially 
only gas on near-plunging high eccentricity orbits arrives in the 
innermost region, but with time gas circularises to e <C 1. 



Particles with Tc < r-acc are in general accreted by Sgr A* in 
our formalism (unless their orbit is changed by interactions 
before they make it inside race). Particles with Vc > race set- 
tle into a thin rotationally supported disc, or into sheared 
out filaments (the remains of the clouds). Note that due to 
the complicated geometry, some gas particles pass through 
shocks several times before reaching their final destination 
(i.e. the disc). 

Figure |3] shows the profiles of orbital eccentricity (de- 
fined on radial shells) plotted for several snapshots. These 
profiles again show that soon after the collision the inner 
part of the computational domain is dominated by gas on 
plunging - high eccentricity - orbits. At later times finite 
net angular momentum of gas and shocks force the gas to 
circularise outside the inner boundary of the simulation do- 
main. It is notable that gas circularises faster at smaller 
radii, which is naturally expected as the orbital time is 
shortest there. Therefore, towards later times, nearly cir- 
cular gaseous orbits are established in the innermost few 
arcseconds, whereas eccentric ones dominate at larger radii. 



4-2.2 Structures in angular momentum space 

The various structures that are formed in the simulation can 
also be identified by the orientation of the angular momen- 
tum vector L of a particle's orbit. We define two angles, 6 
and (j), so that — |L| cos 9, where Lz is the z-component of 
L, and tan(^ = Ly/Lx. We then plot the values of these an- 
gles for individual SPH particles and stars in the Aitoff pro- 
jection. As the number of SPH particles is too large to show 
each particle, we show only a small fraction of these, ran- 
domly chosen from the total number. This procedure keeps 
the resulting maps reasonably clear while also preserving the 
shape of the distribution. 

The result is shown in Figure |4] for a selection of times. 
At t = all the gas particles are concentrated in the two 
clouds, with the finite spread in the distribution arising as a 
result of all the particles within a cloud having the same ve- 
locity whilst possessing a finite spread in positions about the 
centre. The difference in the orientations of the two clouds 
is almost 120° (cf. Table 1). By t = 100, the clouds have 
collided and the gas particles have experienced 'kicks' in 
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Figure 4. Angular momentum direction in terms of the azimuthal, 9, and the polar, <f>, angles. The frames are shown at t = 
0, 100,250, 1000. Gas particles are colour-coded; R < 1" (red), R < 5" (orange), R > 5" (blue). 



diflterent directions, filling the entire range in 6. The (j> val- 
ues are restricted to a rather narrow band running about 
(j) = -1-90° and (p = —90° lines. Note that this band struc- 
ture forms a plane L^^ ~ in L-space. This is expected since 
gas particles experience the collision close to the a;-axis, i.e., 
where y/x ~ z/x ~ 0. Therefore, the a;-component of the 
resulting angular momentum, Lx ~ yVz — zVy is very small. 

The plots are colour-coded to show the behaviour of 
the gas in the inner arcsecond (red), the inner 5 arcseconds 
(orange), and outside the inner 5 arcseconds (blue). Whilst 
the distribution just after the collision is fairly spread out 
for all the three regions, we note a clear difference at later 
times. In particular, the innermost disc occupies a region 
positioned in between the two clouds' original positions. It 
is also the least spread out structure, defining a thin disc 
which is only slightly warped by t = 1000. The innermost 
disc orientation does evolve with time, however, as we noted 
already in i )4.1l since matter infall on the disc continues 
throughout the simulation. 

The gas coloured in orange demonstrates a greater ex- 
tent of warping, and the region outside this (> 5") cannot 
even be classified into a single structure (cf. Figure [!}. Nev- 
ertheless, it is clear that the inner and outer gas distributions 
are similarly oriented with respect to each other, with the 
majority of the outer gas angular momentum distributed in 
6 and (j> between the initial values of the clouds. 



4-2.3 Star formation 

Stars in SI form in both the disc and the primary filament. 
In the disc, stars form first in the inner arcsecond at i ~ 
900 and later at a radius of ~ 5 - 8" at t ~ 1700. In the 
filament, star formation is approximately co-eval with the 
inner arcsecond population and occurs at a radius of ~ 15 — 
25". Mass functions in the disc become top-heavy relatively 
quickly (as the dynamical time is short), with high-mass 
stars (M ~ 100 M0) appearing by t « 1100, whilst stars 
in the filament remain entirely low-mass (~ 0.1 — 1 M©) for 
the duration of the simulation. The disc and filament stellar 
orbits are inclined to each other by ~ 60° , similar to the gas 
orbital inclination. 

All three stellar populations can be seen in Figure [5] 
the "inner disc" population (within the inner 1 — 2"), the 
"mid-range disc" population (5 — 8") and the "filament" 
population. The orbits in the inner disc population are al- 
most circular, e ~ 0.05, whilst stars in the mid-range disc 
have eccentricities of e ~ 0.2. Both of these are therefore 
in goo d agreement wit h the clockwise feature of the obser- 
vations |Lu^^et^ir|(200a) , the orbits of which are only mildly 
eccentric. The orbit of the filament population is also ec- 
centric, with e ^ 0.2. This is in reasonable agreement with 
the corresponding counter-clockwise feature of the observa- 
tions, although the eccentricity of the latter is considered to 
be somewhat higher with e ~ 0.8 jPaumard et al.. 20061 ). 
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4.3 SIMULATION S2 

4-3.1 Gas dynamics 

The simulation S2 is set up identically to SI except for the 
initial velocity of the secondary cloud (cf. Table 1). As the 
result, the angle between the orbital plane of the two clouds 
is greater, 9 ~ 150°, and the collision between the clouds is 
more grazing (i.e. the impact parameter b is greater). 

Smaller fractions of the clouds' mass are therefore in- 
volved in a bodily collision, and as a result the two clouds 
survive the initial collision better. Indeed, comparing the SI 
and S2 simulations in Figure O one sees that in the lat- 
ter simulation most of the secondary cloud actually does 
survive the collision. This cloud forms a clockwise filament 
which was largely absent in simulation SI. Both S2 filaments 
contain a greater amount of mass than their SI counterpart. 

On the other hand, the parts of the clouds involved in 
the direct collision experience a stronger shock in S2 than 
they do in SI as the velocity vectors in the former are almost 
directly opposing. Thus, despite the larger impact parame- 
ter, the stronger angular momentum cancellation in S2 cre- 
ates an inner r < 1 disc as massive as that in SI. However, 
the "mid-range" disc largely fails to form in S2. 

Another difference between SI and S2 is the evolution 
and the structure of the inner gaseous disc. In both simu- 
lations the angular momentum of the gas infalling onto the 
discs evolves with time, and hence the mid-planes of both 
discs undergo a change of orientation in both 6 and (j) (we 
refer to this as 'midplane rotation'). However the infalling 
gas seems to be more intermittent in S2 compared to SI, 
possessing a larger range of angular momentum. Whilst in 
SI the accreting gas largely adds to the existing disc, creat- 
ing one large, relatively warped disc, in S2 the infalling gas 
causes enhanced disc midplane rotation, and at later times 
{t ~ 400), actually creates a second disc around the first 
at a different orientation. Over time, the first (inner) disc 
aligns somewhat with the second one. This process repeats 
at larger radii as gas continues to fall in with a different 
orientation of angular momenta to that of the existing disc. 
The eventual result is that the S2 disc is significantly more 
warped than its SI counterpart. Figure [7] compares the inner 
discs of SI and S2 in the process of forming. 

The difference in gas distribution between the two sim- 
ulations can be also noted by looking at the respective his- 
togram plots of circularisation radius (Figure [2] and Figure 
[8] respectively) . A tail to small rc appears faster in S2, by 
t = 25 rather than f = 50 in SI. The radial profile of gas 
in S2 by t = 500 has a gap between the inner and outer 
distributions, and is far more peaked than the flatter, more 
even distribution of SI at an equivalent time. 

4-3.2 Star formation 

Stars form earlier in S2 than in SI, and they form in the disc 
flrst, by t = 400. The preferential sites of star formation in 
the disc appear to be dense concentric ring structures - see 
Figure [7] - we discuss these in more detail in H4.6I Owing 
to the significantly peaked distribution of gas (Figure |8} in 
S2, there is a far greater amount of star formation at both 
small and large radii compared to SI at an equivalent time. 
This can be seen in Figure [5] The stars formed in the disc 
are again top-heavy and in the filaments are again low-mass. 
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Figure 8. The distribution of SPH particles over circularisation 
radius for simulation S2. The black, cyan, blue and red curves 
correspond to times t = 0, 25, 50, 500 respectively. 

The stars in the filaments show a tendency to form inside 
clusters - this effect is discussed in more detail in M-4\ 

4-3-3 Structures in angular momentum space 

Figure |S] shows the orientation of the angular momentum 
vectors for orbits of gas (in grey) and stars (in colour) in 
the Aitoff projection, for SI and S2 at the same times as 
in Figure [5] We identify stars of different masses using dif- 
ferent colours; green (0.1 — IM©), cyan (1 — 10 M0), blue 
(> 10 M0), and magenta (> 150 M©). As already noted, 
in S2 the secondary cloud largely survives the collision and 
forms a stream, which later collapses into stars. These stars 
can be found at the original orientation of that cloud in the 
right panel in Figure |6] (green at around 9 = —60°). The 
primary and secondary outer stellar populations are there- 
fore inclined at a large angle of ~ 150° to each other, as 
is expected since they are the stellar remnants of the origi- 
nal clouds. The angular momentum vector of the disc stellar 
population is oriented at an angle in between that of the two 
outer stellar populations. 

The other interesting feature of these plots is that the 
stellar disc is strongly warped. The origin of this warp is 
in the warping of the gaseous disc rather than any secular 
evolution once the stars have formed. 

4.4 SIMULATIONS S3 & S4 

S3 and S4 are identical to SI and S2 respectively, except for 
the value of the cooling parameter l3 which is set to a lower 
value of 0.3. 

4-4-1 Gas dynamics 

To show the effect of the faster cooling. Figure |5] contrasts 
gas surface density plots of simulations SI and S3 at time 
t — 150. The bulk dynamics of S3 are the same as in SI, 
but apart from the innermost few arcseconds, maximum gas 
densities are higher in S3 than they are in SI. This is natural 
as faster cooling gas can be compressed to higher densities 
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Figure 5. Projected gas densities and stellar positions for simulations SI (left) and S2 (right) at t = 1955. Stars are shown by coloured 
symbols as following: green Af» = 0.1 — IMq, cyan Af» = 1 — IOMq, blue A/* > 10 Mq, and magenta Af* > 150 M©. Note the survival 
of the secondary filament in S2. The line of sight is along the 2-direction. 




Figure 6. Angular momentum orientations for SPH particles (grey) and stars (colours as in Fig [Sjl in simulations SI (left) and S2 
(right) at t = 1955. Note that the inner stellar disc feature is rather small in angular extent i.e. thin and planar in the left panel, but is 
geometrically thick and significantly warped in the right. 



through shocks and gravity. For the same reason, the fila- 
ments are much better defined in S3 than they are in SI. 
Similar conclusions can be drawn from comparison of runs 
S2 and S4. 



Star formation 

Star formation occurs earlier in S3 than it does in SI. 
In particular, the first stars appear at t = 150 in S3 
and only at t ~ 1000 in SI. This is expected since 
faster cooling facilitates faster gravitational collapse (e.g., 



iGammie. 200]|Rice et al.. 2005h . The S2 to S4 comparison 
is less extreme but still significant, with stars forming in S2 
by t = 400 and in S4 by t = 100. In both (5 = 0.3 runs 
the initial sites of star formation were the orbiting gas fil- 
aments, with the disc stellar populations forming slightly 
later, at t = 350 (S3) and t = 180 (S4). 

The clearly defined concentric ring structures of the SI 
and S2 discs were also seen in S3 and S4, although they did 
not provide the sites of star formation in the latter cases, as 
fragmentation into sink particles occurred some time before 
the rings began forming. In terms of the stellar mass spectra 
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Figure 7. Projected gas density for the inner two arcseconds, showing stars forming (colours as in FiglSj. Left: Simulation SI at 
t = 1000; Right: Simulation S2 at t = 400. The times correspond to when first stars appear in these runs. 




Figure 9. Comparison of projected gas density for simulation SI (left panel, /3 = 1) and S3 (right panel, /3 = 0.3) at t = 150. Note that 
filamentary structures are denser and better defined in S3 than they are in SI. 
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Figure 10. Projected gas density for simulation S4 (top) angular 
momentum of individual gas particles and stars (bottom) at time 
t = 360. In the angular momentum plot, gas is in grey whilst 
stars are colour-coded by mass, as they are in the gas density 
plot, as perO The extended stellar structures in angular momen- 
tum space, such as the linear features seen, correspond here to 
extremely compact "star clusters". In the case of the most ex- 
tended linear feature in the right-hand plot, velocity dispersion 
inside the cluster is comparable to its orbital velocity. 

all the stars formed in S3 and S4 are low mass (~ 0.1 — 
IM0), both in the filament populations and in the disc. 
Again, this is a natural consequence of a short cooling time 
iNavakshin. 200eShlosman fc Begelman. 1989. ). 

4.4-3 Clustered star formation 

The top panel of Figure [10] shows the surface density of 
S4 some time after stars have begun to form. The bottom 
panel shows the orientation of the angular momentum of 
gas and star particles for the same snapshot. In contrast 



to the equivalent plots for SI and S2, one notices extended 
linear structures, "filaments" or "arms" in Figure [TOl These 
structures, extended in angular momentum space, are in fact 
produced by the most compact features in coordinate space 
- the star clusters. These clusters axe seen visually in the top 
panel of Figure 1101 The velocity dispersion of the stars in 
these very dense clusters is high; some are comparable to the 
orbital velocities of the clusters themselves. For example, in 
the bottom panel of Figure IIOI the two longest continuous 
"filaments", separated by 180° in (j), are produced by the 
densest and most massive cluster located at 2: ~ 0, a; = 4, 
y « — 28 in the top panel. Due to its location close to the y- 
axis, the j/-component of the angular momentum of any star 
that belongs to the cluster, ly = —xv^ + zvx ~ 0. Therefore 
the cluster's stars form an incomplete plane ~ in the 
angular momentum space, producing the two (actually one) 
very long "filaments" in the bottom panel of Figure [10] This 
cluster is a remnant of cloud 1, and its mass is ~ 3000 solar 
masses. 

In fact most of the stars in simulations S3 and S4 be- 
long to a cluster. This clustered mode of star formation in 
the fast-cooling runs is due to the gas collapsing promptly 
and thus forming very massive dense gas halos. We sus- 
pect that these results would have been difi'erent had we 
been able to treat radiative feedback from star formation 
in our simulations. We note that in such a dense environ- 
ment, any feedback from star formation activity, such as 
radiation or outfiows, would not have been able to escape 
easily and would therefore have heated the surrounding gas, 
suppr essing further fragmentation and increasing t he Jeans 
mass fNavakshin. 200e' Krumholz fc Bonnell. 200?! ). Rather 
than the multiple low-mass stars in a cluster we would per- 
haps have seen fewer but higher mass stars. The mass spec- 
trum of stars would probably also change if the colliding 
clouds had net rotation (spin) or a turbulent structure be- 
fore the collision, providing some stability against collapse. 
All these efi'ects are however minor as far as orbital motion 
of the stars is concerned, which is the focus of our paper. 

4.4-4 Inner fluctuating disc 

Simulation S4 is the fast-cooling equivalent of S2. As such, 
the inner gaseous disc in S4 is subject to a similar inter- 
mittent infall of gas from larger radii. The newly arriving 
gas tweaks the disc orientation significantly, causing it to 
undergo midplane rotation (refer back to H4.3.1|l . The dif- 
ference here however is that disc is fragmenting into stars 
whilst its orientation is changing. Interestingly, these stars 
do not follow the evolving disc orientation but rather re- 
main in their original configuration. This is clearly seen in 
Figure 1111 where the stellar disc remembers the "old" ori- 
entation of the gaseous disc in which the stars were born, 
whereas the gaseous disc evolves to quite a different orienta- 
tio n. The possibility of this effec t taking place was suggested 
bv lNavakshin fc Quadra (20oi ). These authors found that 
stars will not follow disc midplane changes if these occur 
faster than the "critical rotation time", which is estimated 
to be about 500 code units for the disc in Figure 1111 Our 
simulations are thus consistent with these predictions as the 
disc orientation changed on a timescale of just 80 code units. 

As a result of a several star forming events in gaseous 
discs of different orientations, the inner stellar disc in S4 is 
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Figure 12. Projected gas density for S5 at i = 2696, showing 
stars as well (colours as in Fig[5ll. Note that the stellar populations 
all possess similar values for eccentricities. 

very much diflterent from that found in S2, and even more 
so when compared to SI. In the latter cases, when coohng 
is more gradual (i.e. /3 = 1), stars only form in the inner 
gaseous disc once it has settled into a relatively stable ori- 
entation. Hence the resulting distribution of stars is corre- 
spondingly thin (although in S2 it is significantly thicker 
than SI, cf. Figure |6]). Such a distribution is broadly con- 
sistent with the observed orbit s of young massive st ars in 
the clockwise disc in the GC (|Paumard et al.. 20061 ). Due 
to gaseous disc midplane changes coupled with quicker frag- 
mentation, then, the stellar disc in S4 is much thicker, with 
H/R ~ 1. This particular simulation thus fails to account 
for the most prominent feature of the observational data; 
namely, a thin inner stellar disc. 

4.5 SIMULATIONS S5 & S6 

Simulation S5 is somewhat distinct from all the rest due to 
its small impact parameter, b, between the gas clouds (see 
Table 1). As a result of this fact and the relatively large 
value of /3, a higher degree of mixing is achieved. Figure [121 
shows a snapshot from this simulation, showing both gas 
and stars. Compared with all the other runs, the accretion 
disc extends out to as much as r ~ 20 and represents a very 
coherent disc with little warping. Outer gas filaments do not 
form due to an almost complete mixing of the original gas 
clouds. 

This simulation is an excellent example of how just a 
small change in initial conditions can yield drastically dif- 
ferent results. Comparing the parameters for S5 with those 
for S2 (refer to Table [1} , we see that there are only rela- 



tively minor differences (the largest one being the different 
in impact parameter). However, towards later times the be- 
haviour of the gas diverges considerably, with S5 ending up 
with a single, extended disc and similarly-oriented stellar 
populations and S2 possessing a small, highly warped stel- 
lar disc together with outer stellar populations inclined at a 
large angle from each other. 

Simulation S6 is again different in that it possesses the 
largest collision velocity of all the simulations, whilst retain- 
ing a cooling parameter of (3 — 1. One therefore expects the 
resulting thermal expansion of the clouds to significantly 
modify the velocities of the particles and encourage mixing. 
However, since the impact parameter is relatively large, this 
simulation actually has more in common with S2 than with 
S5, particularly in the inner parts where the gas forms into 
a double disc structure (refer back to i ]4. 3.1(1 . The gas dy- 
namics are therefore largely equivalent to S2, and as such 
we will not discuss this simulation in detail. 

4.6 Star formation in rings 

One particular interesting feature of our simulations is that 
the disc surface density exhibits slightly eccentric "rings" 
in the inner 0.5 — 1". Star formation occurs preferentially 
in these rings (see Figure [7| in the slower cooling runs. 
The inner discs in runs S1-S6 are eccentric, warped, with 
a radial gradient in both of these properties, and evolving 
under the influence of a continuous but variable infall of 
material from larger radii, in a fixed (BH + stellar cusp) 
potential. It is difficult to disentangle these various effects 
from one another to understand the origin of the rings in 
detail, but we believe that they are partly a consequence of 
a fast cooling time (since tcooi oc (JR^^^) coupled with most 
likely a number of the above properties of the gas in the 
inner parts. We base this conclusion on our unsuccessful at- 
tempts to reproduce the rings in cleaner purpose-designed 
simulations. In particular we ran simulations with eccentric 
and circular, fia t or warped discs, but f ound no ring struc- 
tures. Similarly, iNavakshin et al. (2007h did not encounter 
any such structures in their simulations of circular or ec- 
centric self-gravitating discs. We plan to continue to inves- 
tigate the cause of the rings, and in particular to determine 
whether they are a numerical or physical artifact. It should 
be noted at this point that the effect of the rings on our 
model and overall conclusions is (fortunately) negligible, as 
they only affect star formation within the inner ~ 0.5", a 
region in which fragmentatio n would be suppressed through 
accretion luminosity ((Navakshin. 200^ . if we were to in- 
clude this in our model. 



5 RADIAL DISTRIBUTION OF STARS 

Interestingly, the most robust result from all of the 6 runs 
completed is the fact that distribution of stellar mass in an 
annulus, ciM* = Ti,{R)2-nRdR, versus radius roughly follows 
the law 

S*(i?) « ^ . (10) 

We show these distributions for runs SI, S2, S3 and S5 (all 
the simulations that have progressed far enough to trans- 
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Figure 11. The inner 1.5 arcseconds of the simulation S4 at t = 180 (left panel) and t = 260 (right panel). Note that stars born in the 
gaseous disc at time earlier than t = 180 kept their orbital orientation whereas the gaseous disc has evolved due to deposition of new 
matter from larger radii. 



form a significant fraction of the gas into stars). To facili- 
tate comparison with observational data, we nse projected 
radius instead of the proper 3D radius, which is not di- 
rectly known in obs ervations of young stars in the GC (e.g., 
iGenzel et al.. 20031 ). As the viewing angle of the stellar sys- 
tem modelled here is arbitrary, we chose to plot the dis- 
tributions along the three axes of the simulations, so that 
r = \/x'^ -\-y'^ for the black symbols, r — \J'\p' -\- for the 
red symbols, and r = \J for the blue symbols. The 
lines in the Figure show r~'^ (dotted) and r~^l^ (dashed) 
power laws. The latter corresponds to that predicted by non 
self-gravttating standard accretion disc theory. 

This result is remarkable given very large differences in 
the 3D arrangments of the stellar structures in our simula- 
tions. The 1/R^ law is in go od agreement with th e observed 
distribution of "disc" stars l|Paumard et al.. 20061) . We have 
not been able to attribute a simple anaytical explanation 
to this scaling, other than to suggest that cancellations of 
angular momenta in shocks appears to be very efficient in 
driving gas to smaller radii. On the referee's request how- 
ever we have looked at the distribution of the gas for a given 
N-body velocity profile emerging from the initial cloud colli- 
sion, and have attempted to fit the resulting circularisation 
radius distribution to our result for S«(7?). The thermal 
expansion of our clouds as they collide can be modelled to 
a first approximation as an isotropic, delta-function veloc- 
ity profile, which for small values of Vc yields a slope of 
E oc . Introducing a gaussian velocity profile into the 
isotropic distribution, centered about low velocities, acts to 
steepen this slope closer to the desired R~^. Beyond this 
simple N-body model, shocks and small-scale collisions else- 
where in the simulation will drive gas deeper into the po- 



tential well, increasing the slope even further. It would be 
interesting to explore a larger parameter space to quantify 
the robustness of these results. 



6 SGR A* FEEDING DURING AND AFTER 
THE CLOUD COLLISION 

Figure [14] shows the mass accretion rate for the central black 
hole for simulations S1-S4. Formally, the Figure indicates 
that our simulations provide a sustained super-Eddington 
accretion rate (the Eddington accretion rate for Sgr A* be- 
ing ~ 0.03 Mq yr~^), in some cases for between 10'* to 10^ 
yrs. However it must be remembered that we do not resolve 
gas dynamics inside the accretion radius. We expect that 
material will form a disc there, and accretion will proceed 
viscously. The viscous time scale, tvisc, depends on the tem- 
pe rature in the disc midplane and the viscosity parameter, 
a (|Shakura fc Sunvaev. 1973 ). The midplane temperature 
in the inner arcsecond is aro und 10^ K in both the stan - 
dard (non self-gravitating, e.g., lNavakshin fc Cuadra. 20051 ) 
and the self-gravitating regimes ( Navakshin. 20od) . yielding 
H/R ~ 0.01. Hence the viscous time is 



tvisc = 6 X 10^ years Oq J 



R 



WOH 



R 



3/2 



(11) 



where ao.i ~ a/0.1 and R is the radius of the disc in arc- 
seconds. With these fiducial numbers, the viscous timescale 
coincides with the age of the young massive sta rs in the 
GC (e.g.. lKrabbe et al.. 199EtPaumard et al.. 20061 ). As a is 
highly uncertain, we shall then consider the two opposite 
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Figure 13. Stellar surface density, E», versus projected radius, for runs SI (top left, at t = 1955), S2 (top right, at t = 2243), S3 
(bottom left, at t = 1396) and S5 (bottom right, at t = 2696). Different colours show the three different orientations, along z, x and y 
axes (black, red, and blue respectively). Power laws are overlaid for (dotted) and r~^^^ (dashed). All four plots correspond to the 
maximum time to which each simulation has been run. 



cases. If a = 0.1, tvisc ~ a few million years (« 10^ years at 
0.2"), and we expect that gas would have mainly accreted 
onto Sgr A* by now. This accretion rate would be a signifi- 
cant fraction of the Eddington accretion rate. Standard disc 
accretion would then generate as much as 2 x 10^^ erg Macc,3 
of radiative energy, where Mace, 3 is the gas mass accreted by 
Sgr A* in units of 10^ Mq . A similar amount of energy could 
have been released as energetic outflows. There is currently 
no evidence for such a bright and relatively recent period of 
accretion activity for Sgr A*. 

In the opposite limit, i.e. if a ^ 0.1, and tvisc ~ a 
few million years, a gaseous disc should still be present as 
self-consistent modelling pre dicts that stars should not b e 
forming inside i? < 0.3" - 1" (|Navakshin. 200dLevin. 20061 ). 
There are however very strong observational con- 
straints on the absense of such a massive gaseous 



disc in the inner arcsecond of our Galaxy ( e.g., 
iFalcke fc Melia. lOQlNaravan. 200jCuadra et al.. 20031 ). 

Therefore, our simulations seem to over-predict the 
amount of gaseous material deposited in the inner ~ 0.3". 
We take this as an indication that a better model would per- 
haps involve a gaseous cloud of a larger geometric size, thus 
shifting all spatial scales outwards (see further discussion of 
this in W^ . 



7 DISCUSSION 

7.1 Summary of main results 

In terms of the final distribution of gas and resulting stel- 
lar populations, our results can be divided into two main 
categories. 
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Figure 14. The time evolution of the accreted mass onto the central black hole (solid line) and the corresponding accretion rate (dotted 
line) compared to the Eddington limit for Sgr A* (dotted red line), for simulations SI (top left), S2 (top right), S3 (bottom left), and S4 
(bottom right) 



(i) For those runs with a relatively small impact param- 
eter, namely SI, S3 and S5, the initially small-scale disc 
around the black hole grows by steady accumulation of gas 
to extend out to 7? ~ 15". As a result, more than one stellar 
population is seen in the disc. In addition to the ubiquitous 
stars in the inner few arcseconds, a "mid-range" popula- 
tion is seen, at J? ~ 5 — 8". Populations in the filaments 
vary; in SI, only the primary filament forms a stellar pop- 
ulation, as the secondary has accreted onto the disc by the 
time stars begin to form; in S3, although the gas dynamics 
are the same, the secondary filament forms stars before it 
can accrete, resulting in two filamentary populations; and 
in S5, due to the small impact parameter a single large disc 
is created, resulting in three disc stellar populations all at a 
similar orientation. 

(ii) For the runs with a larger impact parameter, namely 
S2, S4 and S6, the small-scale disc grows very little over 
the course of the simulation, staying within a radius of a 
few arcseconds. The feeding of the disc is far from steady, 
generating enhanced midplane rotation, which, coupled with 
star formation (particularly in the case of faster cooling, e.g., 
S4) results in geometrically thick stellar populations. Both 
the primary and secondary filaments survive, and each form 
stellar populations inclined at a large angle to each other. 



The other results of our simulations, presented roughly 
in the order of their commonality to the 6 simulations, are 
the following. Firstly, the formation of a gaseous nearly cir- 
cular disc in the inner region of the computational domain is 
common to all runs, as is the ensuing formation of stars on 
similar near-circular orbits. This is natural as the dynam- 
ical time in the innermost disc is only ~ 60r^^^ years. On 
the scale of the inner disc, r ~ 1, the disc makes tens to 
hundreds of revolutions during the simulation, allowing for 
a near circularisation of the gaseous orbits. Conversely, the 
outer gaseous stream becomes self-gravitating much faster 
than it could circularise (or even become an eccentric disc), 
and hence orbits of stars in that region are more eccen- 
tric. This division on inner near-circular and outer eccentric 
orbits is in br oad agreement with t he properties of the ob- 
served orbits (iPaumard et al.. 2006h . One discrepancy how- 
ever is the eccentricity of the counter-clockwise feature of 
the observations compared to the outer stellar populations 
as seen in our model; the former is believed to possess a value 
of e ~ 0.6 — 0.8 whilst the latter we find is only e ~ 0.2 — 0.4. 
We believe that simulations that use a strongly eccentric ini- 
tial orbit for one of the clouds should be able to match the 
observations more closely. 
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Another robust result of our simulations is that all of the 
gaseous and ensuing stellar discs are significantly warped, 
by between 30° to 60° measured from the inner to the outer 
radii. Interestingly, most recent analysis of the observations 
reveals strong warps in the clockwise disc system (Bartko et 
al 08, private communication). 

We also found that the radial distribution of stellar 
mass closely follows th e observed S« oc 1/R^ profile of the 
disc stellar populations (jPaumard et al.. 20061 ). This was ob- 
served for all the simulations, although we expect these re- 
sults would change if the colliding clouds moved in simi- 
lar directions, significantly reducing the angular momentum 
cancellation in the shock and the thermal "kick" velocity due 
to the shock. Of course, this latter setup would prevent the 
formation of stellar populations at a large angle of inclina- 
tion with respect to each other (one of the most prominent 
features of the observational data), so it is reasonable that 
we have not tested this possibility in our model. 

For the same initial configurations, simulations with a 
comparatively long cooling time parameter, /3 = 1, lead to 
kinematically less dispersed stellar populations than those 
with faster cooling (i.e., /3 = 0.3). As a result, the longer the 
cooling time, the more closely the resulting stellar system 
can be fit by planar systems in velocity space. The innermost 
stellar disc is t hen reminiscent of the observed clockwise thin 
stellar system (jPaumard et al.. 2006h in terms of stellar or- 
bits. Rapid cooling produces clumpier gas flows that lead 
to significant gaseous disc orientation changes during the 
simulations. Rapidly cooling gas flows naturally form less 
coherent discs; such geometrically thick systems are incom- 
patible with the observations. 

Interestingly, with the chosen initial conditions, faster 
cooling promotes survival of the gaseous streams corre- 
sponding to the orbits of the original clouds. These streams 
fragment and form stars mainly in a clustered mode, al- 
though this might be expected to depend on the details of 
the radiation feedback from young stars, which is not mod- 
elled in this set of simulations. 

Our slower cooling runs produce inner and outer stellar 
systems that are inclined to each other by only about 40°, 
which is much sm aller than the observe d angle of ~ 110° 
between the discs (jPaumard et al.. 20061 ). This is a natural 
consequence of angular momentum conservation in a "vis- 
cous" , well-mixed gas flow - the disc orientation in our model 
will always lie in-between the positions of the initial clouds 
in angular momentum space. 



7.2 Origin of the observed stars: one or two 
clouds? 

None of our runs reproduce the observed properties of stars 
and their orbits in detail. This is not surprising as the pa- 
rameter space for initial conditions is quite large, and we 
have explored only a tiny fraction of it here. We feel, how- 
ever, that the agreement with the observations is encourag- 
ing and that the model is promising. To quantify this dis- 
cussion here, we must note that the o bservational picture 
itsel f is still not entirely cl ear. Whereas iGenzel et al. (20031 ) 
and lPaumard et al. (20061 ) classify the orbital stellar distri- 
bution as two disc s, the statistical signiflcance of the second 
disc is disputed bv lLu et al. (200dLu et al. (20081 ). Although 



we tend to agree with the first opinion, we shall consider 
both of these possibilities without bias here. 

Let us first discuss the origin of the clockwise stellar disc 
in the inner parsec of the GC. The existence of this disc is not 
in question since the w ork of iLevin fc Beloborodov (20031 ) 
and lGenzel et al. (20031 ). Concentrating on this feature only, 
a viable model can be simplified to a self-collision of a large 
molecula r cloud rather tha n collision of two clouds. Recently, 
IWardle fc Yusef-Zadeh (200S ). argued that a viable progen- 
itor for the stellar discs is a large-scale molecular cloud, such 
as the 50 km s~^ GMC actually observed near the GC. In 
this model, the cloud infalls almost radially onto Sgr A*, 
temporarily engulfing the SMBH. While most of the cloud 
nearly passes by, gas within one parsec or so gets captured by 
the SMBH and settles into a bound disc of gas. This model 
may be also useful in explaining not only the observed stel- 
lar disc within the inner parsec but also the CND further 
out at ~ 10 pc. 

Even more recently, iBonnell fc Rice (20081 ) have per- 
formed nume rical simulations of a scenario similar to that 
suggested by IWardle fc Yusef-Zadeh (20081 ). In particular, 
an infalling turbulent molecular cloud was set on a sub- 
parsec impact parameter with respect to the SMBH. While 
their numerical method did not include radiative transfer 
directly (it is of course also not included here), they did use 
an effective equation of state that to some degree mimics 
the radiative effects. It was found that the cloud becomes 
tidally disrupted and forms an eccentric disc around the 
SMBH. This disc later fragmented into stars with a top- 
heavy IMF resulting from inefficient cooling. Therefore, the 
one cloud scenario appears to be attractive both qualita- 
tively and quantitatively as an explanation for the origin of 
the clockwise stellar disc. 

Now, considering the rest of the young massive stars 
that do not belong to the clockwise disc kinematically, we 
note that these account for no less than 50% of the over- 
all stellar population. Whether the majority of these stars 
form the second disc or not, their velocity vectors differ from 
the clockwise disc by up to a few hundred km/sec (which is 
the order of circular velocity at these distances). The ve- 
locity dispersion even in a very compact Giant Molecular 
Cloud, with mass of say 3 x 10'* Mq and size of 1 parsec, 
is only ~ 10 km/sec. Thus the only way to create the ob- 
served kinematically distinct population of stars would be to 
postulate the existence of two or more streams (filaments) 
inside the cloud that pass on opposite sides of Sgr A* and 
do not get completely mixed before forming stars. Given 
our numerical experiments in this paper, this does not seem 
implausible if the cooling time is short /3 ^ 1. What is in- 
teresting in this scenario is that the massive stars of the 
counter-clockwise population would then have to form very 
quickly, i.e., on a dynamical timescale, or else gaseous orbits 
would be mixed. The rotation period scales approximately 
as Trot = 3000(-R/5")^'^^ years, so this is quite fast indeed. 

I t should also be noted that recent N-bo dy simula- 
tions (jCuadra et al.. 200jAlexander et al.. 200?! ) imply that 
it would have been very difficult for the high eccentricities 
and inclinations of the dynamically hotter counter-clockwise 
feature to have been formed from a flat, cold disc via scatter- 
ing processes. A single disc progenitor for both GC stellar 
features is therefore largely ruled out. In contrast, in the 
case of the collision of two clouds as considered here, it is 
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almost too easy to obtain an inner near-circular disc and 
a kinematically diverse stellar population farther out. We 
therefore favor a model where a GMC collided with a pre- 
existing cloud or structure, such as a massive larger-scale 
disc, e.g., similar to the observed CND. 

7.3 Size of the cloud(s) 

The observed well defined, flat, geometrically thin and 
near-circular/mildly eccentric clockwise stellar system 
l|Paumard et al.. 20061 ) is best created via a gentle accu- 
mulation of gas. Several independent major gas deposition 
events lead to a warped disc, and/or mixed systems consist- 
ing of several stellar rings or discs co-existing at the same 
radius. To avoid this happening, the inner disc must be cre- 
ated on time scale longer than the critical rotation time, 
which is estimated at tcr ~ few xlO^ years. 

W hile these results are bas ed on analytical argu- 
ments (jNavakshin fc Cuadra. 2005h and our simplified one- 
parameter /^-cooling model, we note that the results of 
iBonnell fc Rice f2008h corroborate this as their inner stel- 
lar disc appears to b e too eccentric to match the data of 
iPaumard et al. f2006l ). 

Deposition of gas in the inner disc takes place on the 
longest of two timescales: the cooling time tcooi and the col- 
lision time, tcou ~ -Rci/wci, where Rd and Wci are the cloud's 
size and velocity magnitude. In the appendix we estimate 
the realistic cloud cooling time during the collision, and show 
that it is always much shorter than the dynamical time un- 
less magnetic fields are very important. 

We are thus left with the only option to require the 
collision itself be more prolonged than tcr- Estimating the 
velocity of the cloud at Vd ~ 150 km/sec, which is of order 
of circular velocities in the inner Galaxy outside the inner 
parsec, we find tcoU = Rd/vd ~ lO^years -Rci,pc, where Rd,pc 
is the size of the cloud in parsecs. We hence require the cloud 
to be larger than a few parsecs to satisfy tcoii ~ ta- Note that 
this size is not necessarily the original size of the cloud if the 
cloud gets tidally disrupted before it makes the impact. In 
the latter case we can take Rd to be the radial distance 
to the centre of the Galaxy at which the tidal disruption 
took place. Finally, the location of the collision should not 
be too far from the central parsec, or else too much angular 
momentum would have to be lost to deposit a significant 
amount of gas at ~ 0.1 pc. 

Another argument going in the same direction comes 
from a comparison of the radial distribution of gas and stars 
in our simulations w ith the observed stellar distribution 
jPaumard et al.. 20061 ). The former is too compact, i.e., all 
of our simulations deposited too much mass within the inner 
arcsecond. In addition, if that was indeed the case 6 million 
years ago, then Sgr A* would have received a significant 
amount of fuel, enough to become at least a bright AGN. 
Given the long viscous times in the inner arcsecond, Sgr A* 
could actually continue to accrete this fuel now. However, 
it is well known that there is no geometrically thin and 
optically thick dis c inside the inner arcsecond of Sgr A* 
l|Falcke fc MeUa. 19 97Naravan. 2002 Cuadra et al.. 200t ). 
Eliminating the gaseous disc by star formation is not 
an option as there are not enough massive young stars 
observed there. 

Taking all these constraints together, we believe that 



the most realistic scenario would be a GMC of the order of 
a few parsecs in size striking the CND at the distance of a 
few parsecs from Sgr A* . This scenario could perhaps explain 
the origin of the inner edge of the CND at ii ~ 2 pc if the 
refilling time scale is longer than the age of the young stars. 
Alternatively such a cloud could self-collide if the impact 
parameter with respect to Sgr A* is small enough, but the 
cloud needs to be very structured, e.g., essentially consist of 
several smaller clouds or filaments. 



8 CONCLUSIONS 

In this paper we presented several simulations of cloud- 
cloud collisions aimed at reproducing gas flows that could 
have formed (a) gaseous disc(s) in the central parsec of our 
Galaxy, as well as the resulting star formation. We found the 
gas cooling time and the impact parameter of the collision to 
influence the outcome significantly. Nevertheless, there are 
several robust results; (a) the inner near-circular and outer 
eccentric orbital structure of the stars formed there; (b) a 
sharply peaked mass distribution of stars E*(i?) ~ 1/-R^; (c) 
that the gaseous and stellar discs are warped. These results 
are in good accordance with the observations. The break- 
down of the stellar system into one or more components is 
sensitive to the initial conditions and also the cooling pa- 
rameter. 

It appears that a GMC with a size of one to a few par- 
secs, self-colliding on a nearly radial orbit, or striking the 
CND at the distance of a few parsec from Sgr A* could 
explain the known observational data satisfactorily. Future 
observations and modelling might put interesting upper lim- 
its on the timescale over which the massive stars formed in 
the GC, and detail the structure and orbit of such a GMC. 
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APPENDIX A: ON REALISTIC ESTIMATES OF 
THE COOLING TIME DURING THE CLOUD 
COLLISION 

Here we try to estimate the realistic cooling time for the gas 
clouds in the initial collision event. Clouds moving at a few 
hundred km/soc would yield a maximum shock temperature 
of a few million K. Using the optically thin cooling function 
at these temperatures (dominated by Bremsstrahlung and 
metal line cooling) and initial cloud densities, one finds that 
the optically thin cooling time is 3 — 5 orders of magnitude 
shorter than dynamical time at R ^ 25", the latter being 
tdyn = 8 X 10^ yrs. With our chosen parameters for the 
initial conditions, the clouds are Compton-thick, i.e. E ~ 
a few tens of g cm~^. The shock would emit soft X-rays 
through Bremsstrahlung, which would have an absorption 
opacity much higher than that of the Thompson opacity i.e. 
kt = 0.4 cm^ g~^. These X-rays would therefore be com- 
pletely absorbed within the colliding clouds. This would not 
lead to a significant increase in the cooling time, however, 
as the absorbed X-rays would be re-emitted as thcrmalised 
blackbody emission at the effective temperature of the ra- 
diation. The latter is estimated to be around 10^ K. The 
clouds' optical depth to such radiation is of the order of 
a few tons, meaning that the clouds would still be able to 
cool effectively via blackbody emission. The cooling time for 
the cloud in its entirety would therefore remain short (using 
the Stefan- Boltzmann law for radiative emission, tcooi would 
again be ~ 3 orders of magnitude shorter than tdyn). 

However, as was pointed out to us by Yuri Levin, these 
considerations completely neglect magnetic fields that are 
most likely present in the pre-coUision clouds. If the radiative 
cooling time is short, then the shock is essentially isothermal, 
and we can expect the gas to be compressed in the shock 
by a few orders of magnitude for our typical conditions. A 
frozen-in magnetic field would then be amplified by simi- 
lar factors, and magnetic pressure by four to six orders of 
magnitude. The strongly increased magnetic pressure might 
present another form of energy and pressure support against 
gravitational collapse. If magnetic flux is then removed from 
the cloud by magnetic buoyancy instability, this happens on 
time scales longer than dynamical, effectively implying large 
values of /3. 

Summarising, a realistic value of /3 during the initial 
cloud-cloud collision would be very small if magnetic fields 
are not important, or large if strong field amplification oc- 
curs due to the shock. 



